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Abstract 

The goal of the present work is to revisit the cranking formula of the vibrational parameters, especially its well known 
drawbacks. The latter can be summarized as spurious resonances or singularities in the behavior of the mass parameters in 
the limit of unpaired systems. It is found that these problems are simply induced by the presence of two derivatives in the 
formula. In effect, this formula is based on the hypothesis of contributions of excited states due only to two quasiparticles. But 
it turns out that this is not the case for the derivatives. We deduce therefore that the derivatives are not well founded in the 
formula. We propose then simply to suppress these terms from the formula. Although this solution seems to be simplistic, it 
solves definitively all its inherent problems. 
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I. INTRODUCTION 



Collective low lying levels of the nucleus are often deduced numerically from the Interacting Bosons Model (IBM) 1] 
or the Generalized Bohr Hamiltonian (GBH) [2]- [3]. Restricting ourselves to the latter we can say that it is built on the 
basis of seven functions: The collective potential energy of deformation of the nucleus, and for its kinetic-energy part, 
three mass parameters (also called vibrational parameters ) and three moments of inertia. All these functions depend 
on the deformation of the nuclear surface. Usually, the deformation energy can be evaluated in the framework of the 
constrained Hartree-Fock theory (CHF) or by the phenomenological shell correction method. The mass parameters and 
the moments of inertia are often approximated by the cranking formula [IJ-Q or in the self consistent approaches by 
other models Most of the self-consistent formulations are based on the adiabatic time dependent Hartree-Fock- 

Bogolyoubov approximation (ATDHFB) which leads to constrained Hartree-Fock-Bogolioubov (CHFB) calculations 
0-[lO| in which the so-called Thouless-Valatin corrections are neglected. It is to be noted that there are several 
self-consistent formulations for the mass parameters in which always some approximations are made (not always the 
same) . Other types of approaches of the mass parameters use the so-called Generator Coordinate Method combined 
with the Gaussian-Overlap- Approximation (GCM+GOA) Recently new methods have again been developped 

[T^-(l3l|. This leads to a certain confusion and the problem of the evaluation of the mass parameters remains (up to 
now) a controversial question as already noticed in Ref. 13]. 

In this paper we will focus exclusively on the mass parameters, especially on the problems induced by the cranking 
formula, i.e. the "classical" Inglis-Belyaev formula of the vibrational parameters. Indeed, it is well-known that this 
formula leads sometimes to inextricable problems when the pairing correlations are taken into account (by means of 
the BCS model). The transition between normal (A = 0) and superfluid phase (A > Ao ~ 0.3MeV) affects generally 
the magic nuclei near the spherical shape under the changing of the deformation [15j . The problem occurs sometimes 
(not always) exactly in these cases for an unpaired system A ~ 0. In that cases the mass parameters take anomalous 
very large values near a "critical" deformation close to the spherical shape. 

This singular behaviour is well-known and constitutes undoubtedly unphysical effect. It has been early found that 
these problems are due simply to the presence of the derivatives of A (pairing gap) and A (Fermi level) in the formula. 
They have been reported many times Q, [Hj]-[I3] in the litterature, but no solution has been proposed. The authors 
of Ref. and [l5| claim that for sufficiently large pairing gaps A the total mass parameter is essentially given by 
the diagonal part without the derivatives, whereas those of Ref. [l7| affirm that the role of the derivatives is by no 
mean small in the fission process and this leads to contradictory conclusions. Other studies [l8| neglect the derivatives 
without any justification. Some self-consistent calculations met also the same difficulties. For example in Ref. [22j |. 
resonances in mass parameters have also already been noticed. As in the present work they were attributed to the 
derivative of the gap parameter A near the pairing phase transition. In short, up to now the problem remains unclear. 
Curiously, one must point out that contrary to the vibrational parameters, the same formulation (Inglis-Belyaev) for 
the moments of inertia does not exhibit any explicit dependence on A and A (as the I-B formula does for the mass 
parameters) and this explains why the I-B formula for the moments of inertia does not meet such problems. This 
difference appears not so natural and is a part of the motivation of this work. All these problems as well as intensive 
numerical calculations led us to ask ourselves if the presence of these derivatives is well founded. If this is not the 
case, their removal should be justified. In fact, the Inglis-Belyaev formula is based on the fundamental hypothesis on 
contributions of two-quasiparticle states excitations. Rigorously it turns out that the derivatives of A and A do not 
belong to this type of excitations and this must explain their reject from the formula. 

The object of this paper is not so much to tell if this model is good or not or to specify the field of the validity this 
model, etc... This study is simply and wholly devoted to a correction of the Inglis-Belyaev formula in the light of its 
fundamental hypothesis. 

II. HYPOTHESIS OF THE TWO-QUASIPARTICLE EXCITATIONS OR THE CRANKING INGLIS- 
BELYAEV FORMULA. 

A. Without pairing correlations 

The mass (or vibrational) parameters are given by the Inglis formula @ , Q : 



Where \0) , \M) are respectively the ground state and the excited states of the nucleus. The quantities Em,Eo are 
the associated eigenenergies. In the independent-particle model, whenever the state of the nucleus is assumed to be 
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a Slater determinant (built on single-particle states of the nucleons), the ground state|0) will be of course the one 
where all the particles occupy the lowest states. The excited states \M) will be approached by the one particle-one 
hole configurations. In that case, Eq. (TTJ) becomes: 

<*l£lO (l\^\k) 

D ij {0 1 ,.,/3 n } = 2H 2 E —^7 f 1 — (2) 

l>\,k<\ \ e l — e k) 

where {0x, .,/3 n } or in short {/?} is a set of deformation parameters. The single particle states \l) , |fc) and single par- 
ticle energies ei, €k are given by the Schrodinger equation of the independent-particle model [19j |. i.e. H sp \ v) = e v \v), 
where H sp is the single-particle Hamiltonian). At last A is the Fermi level 

Using the properties {u\ d/dp = {u\ [d/dp, H sp ] / (e„ — e M ) and [d/dp, H sp ] — dH sp /dp Eq.@ becomes 

(*I^IO<i|^l*> 

D ij {p l ,.,p n } = 2h 2 E — f ^ (3) 

i>x,k<\ (.ei-ek) 

H sp is the single-particle Hamiltonian and A is the Fermi level. 



B. With pairing correlations, hypothesis of the two-quasiparticle excitations states 

It must be noted that in Eq. ([3]) the denominator ei — vanishes in the case where the Fermi level coincides with 
two or more degenerate levels. This is the major drawback of the formula. It is possible to overcome this difficulty by 
taking into account the pairing correlations. This can be achieved through the BCS approximation by the following 
replacements in Eq. (TTJ): 

i) the ground state \0) by the BCS state \BCS) . 

ii) the excited states \M) by the two-quasiparticle excitations states \v, fi) — a^ott \BCS) (here we consider only the 
even- even nuclei). 

hi) the energy Eo by Ebcs and Em by the energy of the two quasiparticles, i.e., by E v + E^ + Ebcs- The BCS 
state is defined from the "true" vacuum |0) by: \BCS) = Ilk [uk + Vka^a^j |0). 

(4) 

Where {u v ,v^) are the usual amplitudes of probability and 

E, 



v = ^{e v -X) 2 + A* (5) 



is the so-called quasiparticle energy. 

As shown by Belyaev [20l | or as detailed in appendix tha above formula can be written in an other form: 

d d 
(H Im) (Ml oflT W) i i Bv 8v 

D u (ft, „ft) . („„„„ + W »* 'ft + a^J^ (6 ) 

Beside this formula, there is an other more convenient formulation due to Bes modified slightly by the authors of 
Ref. [2[ where the derivatives du v /d(ii, du v /d(3j of Eq.© are explitly performed (see also the details in the appendix 
of the present paper): 

o, {A , . m-EE ("i a -§i w w d -§if m h- ^£§«r«? (7) 

here the most important quantity concerned by the subject of this paper is R\ (once again see formula (|A10[) in 
appendix how this quantity is obtained): 

„„ , | 9H sp dX (e„ - A) <9A 

dpi dpi A dpi 
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The two quantities of the r.h.s of Eq. © and ([7J are in the adopted order, the so-called "non-diagonal" and the 
"diagonal" parts of the mass parameters. The derivatives are contained in the above diagonal term R". In other 
papers, the cranking formula is usually cast under a slightly different form. 

All these formulae (U]), ([5]), (O and others are equivalent. 

The derivatives contained in Eq can be then actually calculated as in the Ref. [l5j . with the help of the 
following formulae. 



d\ acp i + bd 



Pi 



dfc a 2 + b 2 
dA bcp i — adp i 
Wi = a 2 + b 2 



(9) 
(10) 



with 



a = ]TA£- 3 , 6 = £(e„ - A)£7" 3 , (11) 

% = E A w i"> ^~ 3 > = X> - A ) w i") ^ 3 ( 12 ) 

These equations can be easily derived through the well known properties of the implicit functions. In the following 
the expression "the derivatives" means simply the both derivatives given by Eq. and (ITU)) . 

In the simple BCS theory the gap parameters A and the Fermi level A are solved from the following BCS equations 
(fl~3|) and (fl"4|) as soon as the single-particle spectrum {t u } is known. 

2 Np 1 



v=i \J (e„ — A) 2 + A 2 



WnrZ.V I- C " - | (14) 

(iVp is the number of pairs of particles in numerical calculations). Of course, from equations (j 131) and (|14[) the 
deformation dependence of the eigenenergies e„(/3) involves the ones of A and A. 

Formally, the solution of Eq. ([r5|) and (TUT) amounts to express A and A as functions of the set of the energy levels 



C. Paradox of the formula in an umpaired system 



It is well known that the BCS equations have non-trivial solutions only above a critical value of the strength G 
of the pairing interaction. The trivial solution corresponds theoretically to the value A = of an unpaired system. 
In this case, the mass parameters given by © or ([7]) must reduce to the ones of the formula ©, i-e. the cranking 
formula of the independent-particle model. Indeed, when A = it is quite clear that: 

E v = yj (e„ - A) 2 + A 2 ->E V = \e v - A| 

u v , v v — > or 1 therefore in Eq.© {u v v^ + u^v v Y -> or 1 

In accordance with the above assumption (uyV^ + u^v^Y = or 1, we can define v and /i in a such way e„ > A 
and e,j, < A therefore E v + E M = \e v - A| + |e M - A| = e„ - A + A - e M = e v - e M 

so that it is easy to see that the non-diagonal part of the right hand side of Eq.© reduces effectively to Eq. ©, i.e.: 

This implies the important fact that in this limit (A — > 0), the diagonal part (i.e. the second term) of the r.h.s. of 
Eq. must vanish, i.e. in other words: 
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2h 2 J2^R»R» when A 

However in practice in some rare cases of the pairing phase transtion this does not occur because it happens that this 
term diverges near the breakdown of the pairing correlations, i.e., in practice for very small values of A(~ 0) (see 
numerical example in the text below). This constitutes really a contradiction and a paradox in this formula. 
In the quantity R" of Eq ([SJ the diagonal matrix elements (v\ dH sp /d(3i \ v) are finite and relatively small, it is then 
clear that it is the derivatives dA/d(3i and dX/df3i which cause the problem. These features have been checked in 
numerical calculations. In this respect, the formulae (JSJ) and (|10D which give these derivatives are subject to a major 
drawback because their common denominator, i.e. a 2 + b 2 can accidentally cancel. Let us study briefly this situation. 
In effect, this can be easily explained because in unpaired situation we must have A ~ 0, involving a ~ in Eq. (jlll) . 
In addition, b is defined as a sum of postive and negative values depending on whether the terms are below or above 
the Fermi level. Therefore, it could happen accidentally that b ~ in Eq ((TTj) involving serious drawbacks or at least 
numerical instabilities. 



III. QUANTITIES SUCH AS A AND A ARE NOT CONSISTENT WITH THE HYPOTHESIS OF THE 
INGLIS-BELYAEV FORMULA. 



A. Basic hypothesis of the Inglis-Belyaev formula and simplification of the formula 



In the independent-particle approximation the contributions to the mass parameters are simply due to one particle- 
one hole excitations. Thus in the formulae or (J3J the particle-hole excitations are denoted by the single-particle 
states k and I. When the pairing correlations are taken into account, the contributions are supposed due only to 
two-quasiparticle excitations states (u, ~p) ^ v\ in Eq. which gives rise to the first term of Eq.©. The second term 
of this formula is due to the derivatives of the probability amplitudes and has to be interpreted as two quasiparticle 
excitations of the type (y,v). However this is not true for all the terms entering into the product of the quantities 
R^,Rj ■ Let us re- focus onto the formula ([7]) in which we will replace in the second sum the quantity A by its 
equivalent from the identity A = 2u v v u E v . After simplification of the coefficient of R"Rj we obtain: 

3 ■' M ~ 2h \2g v (E v +Erf ( ' ~W ' M> M ~W W) + 2h lf—8ET Ri R i (15) 

The fundamental point is in this way it is clear that all the quantities in Eq. (115[) are associated to quasiparticle states 
v and fi except the derivative of A and A. Quantities such as A and A appearing in {R-,Rj) (seeEq. ©) which are 
deduced from Eq. (fT3)) - (fT4|) are due to all the spectrum, they are clearly not specifically linked to these two particular 
states (otherwise indices v and /i should appear with these quantities). Therefore they cannot be really considered as 
contributions due to two quasiparticle excitation states which is the basic hypothesis of the Inglis-Belyaev formula. 
Therefore they cannot be taken into account. 

With this additional assumption, the element R" must reduce to nothing but a simple matrix element: 



Consequently this contributes to simplify greatly the formula (|15[) which becomes: 

D ij {Pu ., /?„} = 2fi 2 EE + ^- — P + 2h 2 J2 ^ v v v f °PJ (17) 

It is to be noted that the missing term (fi = v) in the double sum is precisely the contribution of the simple sum of 
the r.h.s of Eq. (|17[) . Therefore, the formula (|17p can be reformulated in a compact form: 



Da ;{/?!,. ,P n } = 2H 2 YJ2 



dH sp dH sp 
{E v + Erf 



(18) 
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In this more "symmetric" form, this formula looks like more naturally to the Inglis-Belyaev formula of the moments 
of inertia which does not contain dependence on the derivatives of A and A. 

The "new" formula corrects the previous paradox because in the case of the phase transition A — > 0, we will have 
in this limit for the diagonal terms v — /i, u v v v + u v v v — for any v and due to the fact that the corresponding 
matrix element (v\ dH sp /d(3i \v) is finite the second term of Eq. (fTTj) tends uniformly toward zero so that Eq. p7|) 
reduces in this case to the equation of the unpaired system without any problem. 



IV. ILLUSTRATION OF THE APPLICATION OF THE INGLIS BELYAEV FORMULA IN THE CASE 
WHERE THE SINGULARITY OCCURS 

This is illustrated in Fig. [T]by the behaviour of the vibrational parameter Bpp(fi, 7 = 0) as a function of the Bom- 
parameter in the case of the magic nuclei 54 6 Xes2- These calculations have been performed for the both formulae 
and (p~7|) . i.e., respectively with and without the derivatives dX/dp and dA/d(3. The resonance (singularity) 
Bp/? ~ TOOOOOO^AfeV^ 1 occurs near the deformation (3 = 0.09 for the formula with the derivatives. This happens 
even if A is very close to 0. Between (3 = and f3 = 0.15 the formula without derivatives gives small (finite) values 
(-B/3/3 ~ 25h 2 MeV~ 1 ). These very small values of the independent particle model are due to the collapse of the 
pairing correlations. In addition, during the phase transition, i.e., for 0.1 < A < 0.2, the vibrational parameters 
increase up to the important value Bpp ~ BQQt^MeV~ l . We have checked that this is due to a pseudo crossing 
levels near the Fermi level. However, in this respect we have futhermore checked carefully that there is absolutely no 
crossing levels near the singularity. Thus the singularity is not a consequence of a crossing levels as it is often claimed 
[T3 |. As said before the explanation comes from the fact that in Eq. © and (ITU1) the denominator simply cancels. 
This demonstrates the weakness of the old formula ([7]) with respect to that proposed in this paper, that is Eq. (|TT)) . 
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FIG. 1: Neutron contribution to the mass parameters Bpp for the magic nucleus 54 6 Xeg2; The calculations are performed 
within the cranking formula including the derivatives and for the same formula without derivating; Note the quasi divergence 
(singularity) of the version with derivatives near the deformation p = 0.1 



G 



V. CONCLUSION 



In some rare but important illustrative cases the application of the Inglis-Belyaev formula to the mass parameters 
reveals incontestable weaknesses in the limit of unpaired systems A — > 0. In effect, this formula leads straightforwardly 
to a major contradiction, that is, not only it does not reduce to the one of the unpaired system in the case A = (which 
is already a contradictory fact) but even gives unphysical (singular) values. It has been reported in the litterature 
that self-consistent calculations meet also the same kind of problems (see text) . After extensive calculations within 
the Inglis-Belyaev formula, we realized that these problems are inherent to a spurious presence of the derivatives of A 
and A in the formula. This led us to "revise" the conception of this formula simply by removing the derivatives which 
are not consistent with the basic hypothesis of the formula, that is to say with two quasiparticle excitation states. 
This is the reason why our proposal cannot be considered as a simple recipe to the limit A = but as a well founded 
rectification of the formula which is thus no more subject to the cited problems and reduces naturally to that of the 
unpaired system in the limit A — > 0. 

Appendix A: The cranking formula with pairing correlations 

We have to calculate the matrix element of the type (n,m\ 8 /8f3i \BCS) which appears in Eq. Q of the text, i.e.: 

D SR tf\-9fc2V- (BCS\d /dPi\u,y.)( V ^\9 /dPj\BCS) 
^ij \P1> •> Pn} — *'* l^ v ,fj, E v +E^ 

keeping in mind however that the differential operator acts not only on the wave functions of the BCS state but also 

on the occupations probabilities Uk , Vk (of the BCS state) which also depend on the deformation parameter /3i we 

have to write. 

8 ( 8 \ ( 8 



88- \ 88- / V 88 , 

l^i \ y% / wave func \ " l / occup.prob 

We must therefore to evaluate successively two types of matrix elements 

1. Calculation of the first type of matrix elements 

For one particle operator we have in second quantization representation: 
8 \ _ , , d 



88 1 = 4% 

U P* / wave func u P l 

Applying this operator on the paired system and using the inverse of the Bogoliubov-Valatin transformation: 
a v = {u v a v + v u a±),a+ = {u v a+ + v v u-) 
We find: 



\BCS)=Y J {v\ir\riatau\BCS) (Al) 

\°Pi J wave func v ^ °P t 

d 8 
= Y,u.a (H -QQ. I/*) (u u a+ + VuOu^u^an + v u a±) \BCS) = Y^ v ,» H go; |M> [u u a+ + v u otp)v^a^ \BCS) because 

a u \BCS) = 1 
Therefore 

|-) \BCS) =^( V | A | M ) [u v atv aa + \BCS) + v va ^v u a+ \BCS) } (A2) 

u Pi J wave func ,, u Pi 



1 / wave func 

8 

We must notice that for the term i/ = /iwe will have the contribution (v\ —— \v) \u v v v a^a^ \BCS) + \BCS)} 

&Pi 

which is a mixing of a two quasparticle-state with a BCS state. Because the state given by Eq. (|A1[) must represent 
only two quasiparticle excitation, we have to exclude the contribution due to the term v = fj, from the sum of this 
equation. This restriction leads to the following formula: 



-Hz) \BCS) = J2 H ■§*. ^ \ BCS ) (A3) 

^ % / wave func 
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It will be noted that the term VvOvV^a^ \BCS) vanishes for v ^= fx in the r.h.s of Eq. (| A2|) . We then calculate then 
the first type of matrix elements: 

h = (n, m| ]T (v\ A u„^a+a+ \BCS) (A4) 

The above form of the formula suggests that the excited states must be of the form \n, m) — a t a j \BCS) = \k,T). 
We obtain then: 

h = (BCS\ aj a k £^ (u\ A | M ) u^a+o* \BCS) 
d 

We use the following usual fermions anticommutation relations: 

{ a k> a l} = { a t ,<*l}=0, { a k' a t}= S kl 

Thus the quantity between brakets of the BCS sate gives: 
(BCS\ a jak a+a± \BCS) = (5^5 vk - &p k 6 ul ) 
We obtain: 

Q Q _ Q 

h = Ytv^fj, W flo" I J") u v y n {^i^uk - &pk8 v j) = (k\ ^5- |0 UfcW; - (Z| — |fc) u 7 U|t with k ^ I Because indexes of 
brakets must be different in Eq (IA4l) . 

Noting that if T is the time-reversal conjugation operator we must have for any operator O 
(p\6\q) = (fp^Of-^fq 

Applying this result for our case and assuming that d/dfti is time-even, i.e. T (d/dfti) T^ 1 — d/dfti, we get: 



I x 'dpi 



Moreover, using the usual phase convention 
Uj = Ul,V£= -v k 
we deduce : 

d d d 

h = (k\ \l) u k vi + (k\ — \l) mv k = (u k v t + uiv k ) (k\ — \l) 

Taking into account that the brakets states in Eq (|A3I) must be different, the final result for I\ given (|A4[) will take 
the following form: 

Il = ( k ~ l \(^-) \BCS) = {u k v l +u l v k ){k\^-\l) with k^l (A5) 

V a Pi J wave func "Pi 

Let be H sp the single-particle Hamiltonian and 



H ' = T,v,n H ( H s P - A ) Im) a+a a - GY^^yQa+a^ajza^ 

the nuclear paired BCS Hamiltonian with the constraint on the particle number. Writing this Hamiltonian in the 
well-known quasiparticles representation H' = Ebcs + ^Z v E v (x^ a. v + residual qp interaction, neglecting (as usual) 
the latter term and using Eq. (|A3|) it is quite easy to establish the following identity 



1 / wave func 



\BCS) = -(k:i\{^£) \BCS) = (E kJ -E BCS )( k ;i\(^ \BCS) 

\ ' / wave junc \ "t>/ wave func 



wave func 



where the eigenenergies E k j corresponding to the excited states fe, l) are given by E k j — Ebcs + Ek + 
so that: 



(k,l\ (^pl \BCS) (k,T\ (^pj \BCS) 

C \ 1 -nsrvsiy. \ f 11 J wave func V < % / wave func 



d P iJ wave func ' ^ E^j - E B CS E k + E t 



8 



Due to the fact that the pairing strength G does not depend on the nuclear deformation, it is clear from the expression 
of H'im the particles representation) that dH'/dfc = d(H sp - A) /dft Therefore (k,7\ (8H' / 'dp t ) wave func \BCS) = 
(k,l\ {dH sp /dpi) \BCS) - d\/dpi{k,l\BCS) = (kj\ (dHsp/dfr) \BCS) 
Here we have (k, I \BCS) = because excited states and bcs state are supposed orthogonal. 

Again using the second quantization formalism (dH sp /dPi) = Ylujtu ( u \ {^H sp / 'dPi) a+a^ and performing then 
exactly the same transformations as before for J]^^ (v\ d/dfii a^a a but this time with ( u \ dH sp /d/3i |/Lt) a+a^ 

we will obtain in the same manner a new form for Eq. (| A5|) : 

h=h{k,t) = (k 1 l\(^) \BCS)=- {Uk l l+U ' Vk) (k\^\l) with k*l (A6) 

V °Pi J wave func h k + OPi 



2. Calculation of the second type of matrix elements 



Recalling that the BCS state is given by: \BCS) = Ilk \Uk + v kO-t a %j |0) ana ^ differentiating this state with respect 
to the probabilty amplitudes, we obtain: 

4g) \BCS) = E T (|J + IL^K + v k aia+) |0) 

^ % / occup.prob ' % ' 1 

We use the evident property: 

rw«* + i°) = ^ + v ^t4r x \bcs) 

Therefore: 



_d_ 

dp* 



\BCS) = £ 7 



occup.prob 

Making an expansion of the inverse operator in a^a^ 

w) lBCS) = E 

" l / occup.prob 

using the inverse of the Bogoliubov-Valatin transformation 



\BCS) 



{-Qj£ + ata£)u T L {1 - v T u T L a+af + (v T u T L a+a±) + ...) 



\BCS) 



7+ - 



(u T a+ + v T a-r) 



We find for the quantity a+ai 

u T a+ + v T a-)(u-a± + u-a T ) 
u T u-a+af + u T vr<^t a r + v T u-a-af + v T v-a-a 



replacing in the above expression and retaining only two creation of quasiparticles with at most products of two 
amplitude probability: 



d 



QR . 

" % / occup.prob 

Noting that: u— 
_d_ 

Oh 



\BCS) = £ T 



\BCS) = Er 



x+a± \BCS) 



—— u T [~v T u T L u T u-oq:c£ f 

r , we find 

dv T du T 

Ut W~ Vt W 

occup.prob L r J% r^i 

The excited states will be necessarily here, of the following form: 
\M)=a+a+\BCS) = \m,m) 
We have therefore to calculate: 

h = {BCS\a w a. m {u m -^- - v m -^)a+a±\BCS) 

due to the normalisation of the excited states, we obtains: 
_ dv m _ du m 

OPi dpi 

knowing that the normalization condition of the probability amplitudes is: 

«m + u m = 1 

we find by differentiation 

dPi dp* 
combining these two relations, we obtain in I^: 

_ 1 du m 

v m dPi 
then, the second term reads: 



-UT-u^al 



\BCS) 



9 



\ ^ / occup.prob U m VPi 

which can be cast as follows: 



(m, m 



h = h{k,l) = (k~l\ 



_9_ 

dpi 



\BCS) =- — ^ wtft k = l 

occup.prob 



(A7) 



The two matrix elements I\ given by Eq. (|A6p and /2 given by Eq. (|A7|) . They correspond respectively to the 
non-diagonal k ^ I and diagonal part k = I of the total contribution. Reassembling the two parts /i and 1% in only 
one formula, we get: 



h+I 2 = (k,l\-^-\BCS) 



(u k vi + uiv k ) dH s 



10(1-4,0 



1 du k 



E k + E t d/3i v k dPi 

Replacing this quantity in Eq. (Q| of section [TT1 noticing that the crossed terms (-Z1/2 and I2I1) cancel in the product 
we find: 



A,{/3i,.,M = 2^ ( " feU ' 



■ uiv k ) 



k.l 



{E k + Ei) 



dH 



Sp 



dp, 



\k) (k\ 



OH, 



I0(i 



1 1 du k du k 
2E k vidl%~dB~ 3 



(A8) 



The expression 



E l 1 du k du k 
2E k v\ dp, d/3i 



(A9) 



meet in the second part of the r.h.s of the above formula (|A8|) can be further clarified. Recalling that the probabilty 
amplitudes are: 

u k = (1/V2) (l + W^TfA^) V2 and v k = (l/V2) (l - e k /^ejT^) ^ 

where:£fe = e k — X is the single-particle energy with respect to the Fermi level A, e k being the single particle energy. 
Since the deformation dependence in u k appears through e k , A, and A, a simple differentiation of u k with respect to 
Pi leads to: 



du h 



1 



-1/2 



dpi 2V2 I y/e 2 k + A 2 



9_£fc 

dpi 



(4 + A 2 



-1/2 



- £fc (e 



1 



multiplying by — and simplifying we get: 

v k 

1 f A de k 9A1 
Wfe 9/3, 2 ( £ 2 + A 2 ) \ dp, £k dp t J 
using e k = e k — A, we obtain explicitly: 



1 9u 



1 



du k 

dpi ~ o 1^2 



1 



v k — 2 (ej* + A 2 ) 
Moreover, noting that 

we find: 

A 



A 



de k 
dp, 



A 



gA 



(c* - A) 



dP, J 



1 ^ u fc 

^9& = 2 (e 2 + A 2 ) 
the quasiparticle energy is E k = (ef 



dH s 

dp. 



dX 

dP 



(e k - A) dA 



1 9u t 



u fe dp., 
where we have put 



A 



9/3, 



9A 

9/3, 



A 9/3, 

Ay so that: 

- A) 9A 1 A 



A 9/3, 



2^ 2 



A 2 



-3/2 



9gfc 
9/3, 



9A 

Wi 



R 



dH sp dX (e fc -A) 9A 

9/3, 1 ; + 9/3, A 9/3, 



(A10) 



Using the result 



10 



1 du k A 



v k dfc " 2E\ 

the product of the similar terms of Eq. (|A9I) gives finally: 
_\_\_duk_du^ _ 1/1 du k \ ( 1 

2^ fe ^ 9/3j - i^k 2Ek y Vk dfii j y Vk df3 . 

~ ^ k 2E k { lX 2El) { 2£f )-^ k 8El UiH i 

The cranking formula of the mass parameters becomes finally: 

D l3 ., M = 2K* ^l^f d\ ^ I*) (*| ^ 10 (1 - + 2^ £ ^fl? (AH) 
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